Realizing Fractional Chern Insulators with Dipolar Spins 
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Strongly correlated quantum systems can exhibit exotic behavior controlled by topology. We 
predict that the v — 1/2 fractional Chern insulator arises naturally in a two-dimensional array of 
driven, dipolar-interacting spins. As a specific implementation, we analyze how to prepare and 
detect synthetic gauge potentials for the rotational excitations of ultra-cold polar molecules trapped 
in a deep optical lattice. While the orbital motion of the molecules is pinned, at finite densities, 
the rotational excitations form a fractional Chern insulator. We present a detailed experimental 
blueprint for 40 K 87 Rb, and demonstrate that the energetics are consistent with near-term capabili- 
ties. Prospects for the realization of such phases in solid-state dipolar systems are discussed as are 
their possible applications. 



The quest to realize novel forms of topological quan- 
tum matter has recently been galvanized by the theoret- 
ical prediction and subsequent experimental observation 
of topological insulators [1, 2]. Such materials harbor an 
insulating bulk, but owing to nontrivial bulk topology, 
they are also characterized by robust conducting surface 
states. Recent theoretical work has shown that combin- 
ing single particle topological bands with strong interac- 
tions can yield so-called fractional Chern insulating (FCI) 
phases [3-8]. Particles injected into these exotic states 
of matter fractionalize into multiple independently prop- 
agating pieces, each of which carries a fraction of the 
original particle's quantum numbers. Unlike traditional 
bosons or fermions, these anyonic excitations accumulate 
a nontrivial phase under exchange. 

While similar effects underpin the fractional quantum 
Hall effect observed in continuum two dimensional elec- 
tron gases [10], fractional Chern insulators, by contrast, 
are lattice dominated. They have an extremely high 
density of correlated particles whose collective excita- 
tions can transform non-trivially under lattice symme- 
tries [8, 9]. In this paper, we predict the existence of a 
FCI state in dipolar interacting spin systems (see Fig. 1). 
This state exhibits fractionalization of the underlying 
spins into quasiparticle pairs with semionic statistics [11]. 
The predicted FCI state may also be viewed as a gapped 
chiral spin liquid [11, 12], a state which has never been 
observed in nature. Such a state cannot be realized in 
conventional electron gases. 

Several recent studies have conjectured the existence 
of fractionalized topological phases in idealized lattice 
models that require sensitively tuned long-range hop- 
ping and interactions [5-7, 13-15]. Broadly speaking, 
two single-particle microscopic ingredients are required, 
both of which find close analogy in the physics of the 
electronic Hall effect. First, just as electrons in a Lan- 



dau level have no dispersion, the dispersion of the lattice 
band-structure must be quenched relative to the energy 
scale of interactions. This enables interactions between 
particles to dominate over the kinetics of their environ- 
ment [13-15]. Second, the flat lattice band should pos- 
sess a non-trivial Chern number, reflecting the underly- 
ing Berry phase accumulated by a particle moving in the 
band-structure. In the context of electronic systems, this 
corresponds to the well-known quantization of the Hall 
conductance cr xy . To this end, the particles in a Chern 
insulator need also be coupled to a background gauge 
potential that breaks time-reversal. Finally, to observe a 
fractionalized insulating state, one must partially fill the 
topological flat band- structure with interacting particles; 
since the FCI state generally competes with superfluid 
and crystalline orders, the resulting phase diagram nat- 
urally exhibits both conventional and topological phases 
(Fig. 1). Up to now, it has been unclear whether such 
exotic fractional Chern insulating phases can be realized 
in any real-world physical system. 

We consider a two-dimensional array of tilted, driven, 
generalized spins interacting exclusively through their in- 
trinsic dipolar interaction, as depicted in Fig. 1A. This 
interaction mediates the long-range hopping of spin-flip 
excitations, and both of the single-particle ingredients 
sketched above arise naturally in this context. The 
quenching of the spin-flip band-structure owes to the 
anisotropy of the dipole-dipole interaction, which yields 
interference between different hopping directions [16]. 
The production of a synthetic background gauge poten- 
tial is accomplished via spatially varying electromagnetic 
radiation [16, 17]. Together, the dipolar anisotropy and 
radiation induce orientation-dependent Aharonov-Bohm 
phases that ultimately generate topologically nontrivial 
flat bands [16]. 

Our approach is applicable to generic pseudospin de- 
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FIG. 1. Realization of a fractional Chern insulator. (A) Schematic representation of the two-dimensional array of polar 
molecules dressed by optical beams (red arrows). Each polar molecules is characterized as an effective pseudo-spin- flip, which 
can hop and interact mediated by the long-range dipolar interaction; <j> represents the Aharanov-Bohm phase which the spin- flip 
acquires as it traverse a plaquette. (inset) Molecules occupy the {A, Y} plane and the rotational quantization axis is set by an 
applied electric field along the z direction. Go and $o define the {x, y, z} axes with respect to the lattice coordinates {A, Y, Z}. 
(B) We consider the J = 0, 1 manifold of each molecule with the |0,0) state representing spin-down. The spin-up state is 
created via optical Raman dressing in the M-conflguration. The optical radiation admits a single dark eigenstate, which is a 
linear combination of the three states in the J — 1 manifold. (C) Phase diagram for 40 K 87 Rb molecules at half-filling with a 
total of N s = 24 sites as a function of electric field strength and tilt Bo- Each phase finds a direct analogy in the language 
of frustrated magnetism and the equivalent nomenclature is given below. The dotted line at \E\ = signifies the fact that a 
minimal electric field is always required to split the degeneracy within the J — 1 manifold. 



grees of freedom, including singlet-triplet states associ- 
ated with solid-state electronic spin dimers, hyperfine 
coupled electronic and nuclear spins and cold atomic sys- 
tems [18-20]. To be specific, here, we focus on an imple- 
mentation using ultra-cold polar molecules trapped in a 
deep two-dimensional optical lattice. Such an implemen- 
tation has many advantages, including local spatial ad- 
dressing, stable long-lived spins, high phase-space density 
and intrinsic dipolar interactions [20, 21]. The molecules 
are subject to a static electric field E tilted with respect 
to the lattice plane (Fig. 1A). We assume that the molec- 
ular motion is pinned, and hence, restrict our attention to 
an effective rotational degree of freedom on each site. In 
particular, we focus on the four lowest rotational levels: 
|0, 0), the rovibrational ground state and the three states 
within the J — 1 manifold (|1,— 1), |1,0), |1,1)), where 
J characterizes the rotational angular momentum of the 
molecules. Here, the quantization axis, £, lies along the 
applied electric field and |J, m) denotes the state adia- 
batically connected (via E) to the rotational eigenstates 
[22]. Each molecule is driven by optical radiation, which 



couples the three J = 1 states to a pair of molecular 
excited states \ei) and |e2), in the so-called M-scheme 
(Fig. IB). The Hamiltonian for each molecule with the 
laser on has the form 



H r = h[\e 1 )(Sl 1 (l, 
+ |e 2 )(fi 3 (l,0|- 



-i| + n 2 (i,o|) 
n 4 <i,i|) + h.c] 



(1) 



in the rotating frame, where Qi are Rabi frequencies serv- 
ing as the control parameters, and h.c. represents the 
hermit ian conjugate terms. The above Hamiltonian ad- 
mits a unique "dark" eigenstate, |t) = i(^2^4|l,— 1) — 
^i^ 4 |l,0) + ^1^3 1 1, 1)), which is decoupled both from 
the excited states and from the radiation field (Q is a 
normalization). Together with the rovibrational ground 
state, which we label as |^), this forms an effective two- 
state spin degree of freedom on each site [16, 22-24]. 

Individual molecules interact with one another via elec- 
tric dipole-dipole interactions, 



di-d i -3(d i -R <i )(d i -Ri i ) , (2) 



3 



where n — l/(47reo) and R^- connects molecules i and 
j (with dipole moment operators and dj). We let 
d be the permanent molecular dipole moment and Rq 
be the nearest-neighbor lattice spacing; by ensuring that 
the characteristic dipolar interaction strength, kcP /Rq, is 
much weaker than the optical dressing, f^, all molecules 
remain within the Hilbert space spanned by , 
Moreover, this interaction is also much weaker than the 
bare rotational splitting 2B (Fig. IB) and thus cannot 
cause transitions that change the total number of |t) exci- 
tations. This effective conservation law suggests the util- 
ity of recasting the system in terms of hardcore bosonic 
operators, a\ = |t) which create spin- flip "parti- 
cles" . Mediated by the dipolar interaction, these molec- 
ular spin-flips hop from site j to site i with amplitude 
tij = — \ Hdd As each hardcore boson harbors 

an electric-field induced dipole moment, there also exist 
long-range density-density interactions of strength = 

(tit? | Hdd + (Uij \ Hdd \Uij) - Hdd ~ 

{\4\j\HddW4\j)' I n combination, this yields a two- 
dimensional model of hardcore lattice bosons, 

H B = ~J2 ^3 a \ a j + V i3 n i n h ( 3 ) 

whose total number, N = Y2i a l a ii 18 conserved [16]. 
Variations in the dipolar- induced on-site potential, ta, 
can be regulated via tensor shifts from the optical lattice 
(see supporting information for details). 

To ensure that our hardcore bosons reside in a topo- 
logical flat band, we adjust the optical beams that dress 
the molecules to produce a square lattice with four types 
of sites, {a, 6, A, B}, as shown in Fig. 2 A (see supporting 
information for details). Owing to interference between 
the dressing lasers, the dark state on each of the sites is a 
different linear combination of the three J = 1 states, im- 
plying that the hardcore boson, aj, is site-dependent. De- 
spite the existence of four unique lattice sites, so long as 
and Vij remain invariant under the direct lattice vec- 
tors g\ and (Fig. 2A), the Hamiltonian retains a two- 
site unit cell. Thus, computing the single-particle band- 
structure produces two bands in momentum space, as 
shown in Fig. 2B. The nonzero Chern number, C = — 1, 
of the bottom band reflects the breaking of time-reversal 
arising from the asymmetry between the intensity of left- 
and right-circularly polarized light fields [16]. 

To characterize the single-particle dispersion, we com- 
pute the flatness ratio, /, between the band-gap and the 
width of the lowest band [13-15]. Numerical optimization 
of the electric-field and the optical dressing yields a vari- 
ety of flat bands with flatness ratio, / > 10. Since an ap- 
plied DC electric field changes the transition dipole mo- 
ments between single-molecule eigenstates, flatness opti- 
mization must be re-performed at each field strength (if 
one desires to retain a flat-band). The optimized band- 
structure depicted in Fig. 2B has / ~ 11.5 and is obtained 
at weak DC electric fields, just strong enough to split the 
degeneracy within the J = 1 manifold (relative to the 




FIG. 2. Topological Flat Bands. (A) Schematic represen- 
tation of the 2D dipolar array, a, b, A and B sites are charac- 
terized by dark eigenstates that are different linear combina- 
tions of the three J — 1 states. Square plaquettes are charac- 
terized by a time-reversal breaking flux, 0, which is staggered 
throughout the lattice. The lattice harbors a two-site unit 
cell and is invariant under translation by direct lattice vectors 
gi and eft. (B) An optimized band-structure in the reduced 
Brillouin zone (RBZ) depicting a flatness ratio / ~ 11.5. The 
lowest band carries Chern index C = — 1. The electric field 
tilt is {Bo, $0} = {0.68, 5.83} and the electromagnetic driving 
parameters are detailed in the supporting information. 



dipolar interaction strength) and to set the quantization 
axis [21]. In addition to the experimental simplicity of 
working with weak DC fields, such a scenario also effec- 
tively eliminates direct long-range interactions between 
the hardcore bosons, as the induced dipole moments of 
|t) and \\) are then negligible. 

With topological flat bands in hand, we now consider 
the actual many-body phases which arise at finite lat- 
tice filling fractions v (number of spin flips per unit cell) . 
To this end, we perform exact diagonalization of the full 
many-body Hamiltonian at v = 1/2 on systems of vary- 
ing sizes up to N s = 44 sites with periodic boundary 
conditions. For weak electric fields tilted near the so- 
called magic angle, 0o = cos _1 (l/v // 3) where the strong 
phase anisotropy between X- and Y- hops produces par- 
ticularly flat Chern bands, diagonalization reveals the 
existence of a bosonic v = 1/2 fractional Chern insu- 
lator. This phase exhibits gapped fractionalized quasi- 
particles of effective charge 1/2 with respect to the mi- 
croscopic spin- flip excitations. As numerical diagnostics, 
these topological features require the presence of two- 
fold ground state degeneracy on a torus (Fig. 3A) and 
a neutral spectral gap that is stable as the system size 
increases (Fig. 3B). The quantity analogous to the Hall 
conductance, a xy = ^ J J F(Q x ,6 y )dQ x dQ y = —0.5, ap- 
pears unambiguously in the response of the system to 
boundary-condition twists {Q x , y } (equivalent to flux in- 
sertion) in the form of a well-quantized many-body Berry 
curvature, F(0 X , B v ) = Im((^|^) - <£|^». 

The counting statistics of low energy quasihole states 
provide a direct diagnostic of the fractionalization of re- 
moved particles. In the continuum, this counting can 
be understood through a generalized Pauli principle, in 
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FIG. 3. Evidence for v = 1/2 FCI state. (A) Exact diagonalization of the full Hamiltonian at v — 1/2 with a total of 
N s = 24 sites and Nb = 6 hardcore bosons. The electric field and driving parameters are identical to those used in Fig. 2b. To 
avoid self-interaction, we truncate the dipolar terms at order l/(3i?o) 3 - There exist two degenerate ground states in momentum 
sectors (k x , k y ) = (0, 0) and (k x , k y ) — (— 7r, 0) consistent with = 1/2 FCI state on a torus (k x , k y are crystal momenta). (B) 
Finite size scaling from N s — 16 to N s — 44 suggests a stable spectral gap in the thermodynamic limit. (C) Quasihole counting 
for the same parameters as in (a) with a single boson removed (Nb = 5). There exists a clear gap below which there are 36 
low-energy quasihole states, consistent with our analytical counting formula [3]. (D) Filling fraction as a function of chemical 
potential, fi — EN b +i — EN b (where EN b is the ground state energy with Nb bosons). Below v — 1/2, there exists evidence of 
a clear compressible superfluid state, while at v — 1/2, there exists a plateau indicative of an incompressible quantum liquid. 
This plateau can also be interpreted as a magnetization plateau in the language of frustrated magnetism [11? ]. (inset) Depicts 
dE/dSo as a function of tilt Go- Phase transitions between the superfluid, Chern insulator and striped solid are evidenced as 
jumps in dE/dOo. Curves from top to bottom are for increasing electric field strength from E = 0.4 — 8 kV/cm. 



which there should be no more than one particle in two 
consecutive Landau-level orbitals; interestingly, the same 
overall counting holds for the lattice-dominated Chern 
insulator [3]. Counting the total number of admissi- 
ble quasihole arrangements on a torus (for a v — 1/2 
FCI) yields, Q torus = fc^) - fc'i^J' where 
N uc = N s /2 is the number of lattice unit cells and 
is the number of hardcore bosons. In our system, nu- 
merically counting the total number of quasihole states 
matches the above formula. For example, the energy 
spectrum of the system (depicted in Fig. 3A) with a single 
boson removed (N uc = 12, = 5) is shown in Fig. 3C; 
from the analytic counting, one obtains Qtorus — 36, in 
precise agreement with the number of states below the 
spectral gap. 

Remaining at v = 1/2, we now probe the many-body 
phases which arise as one varies the DC field strength and 
the tilt, Oo, while adjusting the optical parameters to 
keep the local dark states fixed (see supporting informa- 
tion for details). Changes in the tilt alter the geometry 
of the dipoles and introduce additional dispersion into 



the single-particle bands. On the other hand, increasing 
the electric field strength enhances the long-range inter- 
actions. These qualitative differences in the microscopies 
yield a rich phase diagram exhibiting both conventional 
and topological phases, as shown in Fig. 1C. In addi- 
tion to the FCI phase, there exist four distinct crystalline 
phases at strong DC fields and a large superfluid region 
at moderate fields. While we use the language of lattice 
bosons above, we note that the FCI phase may also be 
interpreted in the language of frustrated magnetism as 
a chiral spin liquid while the competing superfluid and 
crystalline phases correspond to XY ordered magnetic 
and spin density wave phases [11, 25]. As summarized in 
Table L, we characterize the nature of these phases via 
four diagnostics: i) ground- state degeneracy, ii) spectral 
flow under magnetic flux insertion (superfluid response), 
iii) real-space structure factor (n(i?)n(0)), and iv) <j xy , 
the many-body Berry curvature (see supporting informa- 
tion for details). 

Next, we consider a possible route to preparing the v = 
1/2 fractional Chern insulator. In current polar molecule 
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TABLE I. Diagnostics of Many-body Phases 



Phase 


Degeneracy 


SF Response 


Structure Factor 




v — 1/2 Chern insulator 


2 


none 


fluid 


-0.5 


Superfluid 


1 


isotropic 


fluid 


gapless 


Striped supersolid (SSS) 


3 


uni-directional 


stripes 


gapless 


Knight's Move solid (KMS) 


4 


none 


knight's move 





Checkerboard solid (CKB) 


2 


none 


checkerboard 





Striped Solid (SS) 


4 


none 


stripes 






experiments, the spin-flip "vacuum" , corresponding to all 
sites in the |^) state, may be prepared with high fidelity 
from Feshbach molecules by two-photon stimulated Ra- 
man adiabatic passage [21, 26, 27]. Moreover, rotational 
excitations of the molecules can be deterministically in- 
troduced by electromagnetically driving the |^)-to-|t) 
transition, enabling the straightforward preparation of 
finite density, low-temperature states within the super- 
fluid and crystalline phases. A detailed variational mean- 
field study confirms this and generates a non-topological 
phase diagram consistent with the exact diagonalization 
(see supporting information). The topological order of 
the FCI state impedes direct preparation by an analo- 
gous procedure. If the phase boundaries surrounding the 
FCI state are second order, one might attempt to prepare 
this state by adiabatically tuning the electric field across 
the transition. The only known second-order transition 
between a superfluid and the v = 1/2 FCI phase is mul- 
ticritical [28] , which suggests that this phase boundary is 
generically first order, consistent with the numerics pre- 
sented in the inset of Fig. 3D. On the other hand, continu- 
ous Mott insulator to FCI transitions are less finely tuned 
and may constitute a promising avenue for preparation. 
In particular, the striped solid phase may be reduced to 
a simple non-translation-symmetry breaking Mott insu- 
lator in the presence of a one-dimensional superlattice 
potential. This observation is consistent with the exis- 
tence of a weaker, finite-size, cross-over at the FCI to 
striped-solid phase boundary (inset Fig. 3D). 

The most direct experimental realization of our pro- 
posal would be in molecules with a 1 S ground state and 
no hyperfine structure. With quantum degenerate Bose 
gases of multiple nuclear-spin-free isotopes of Sr and Yb 
readily available [29] , such hyperfine- free SrYb molecules 
seem likely to be realized in the near-future. However, 
our proposal can also be carried out in currently available 
ultracold polar molecules, such as 40 K 87 Rb [20], 7 Li 133 Cs 
[27], and 41 K 87 Rb [26]. Here, we focus our discussion 
on 40 K 87 Rb. For the optically excited states |ei) and 
|e2), we propose to use the \ J f ,m f ) = |2, ±2) rotational 
states of the v' = 41 vibrational level of the (3) 1 S + 
electronic state. These states harbor a strong 640 nm 
transition to the ground state [21, 26]. To realize our 
proposal, we require a hierarchy of energy scales corre- 
sponding to, Hiatuce < H hf < fii <C #1,0 - #1,1, where 
Hiattice describes the optical lattice potential and Hhf 



characterizes the molecule's hyperfine structure [21]. For 
40 K 87 Rb, this hierarchy is easily realized since Hhf ~ 1 
MHz, while E li0 - ^1,1 = 160 MHz at a moderate DC 
field strength, E = B/d « 0.5kV/cm. By ensuring that 
the optical dressing (f^) is weak relative to the splitting, 
— we can employ frequency selection during 

the creation of the M-scheme; meanwhile, the condition 
Hiattice % Hhf <C Vti allows us to consider hyperfine 
and tensor-shift effects only after the dark state 
is already defined. Since the composition of the dark 
state is different on various sites, the spin- flip particle 
will be subject to a site-dependent lattice potential and 
hyperfine structure. The application of a magnetic field 
~ 10 3 G (already present in experiments to tune the Fes- 
hbach resonance [20]) defines the decoupled nuclear spin 
basis as an approximate eigenbasis of both (^\Hhf\\) 
and (\\ Hhf \\). It is then possible to find a pair of nu- 
clear spin eigenstates that have good overlap and whose 
compensated energy difference (taking into account Hhf, 
Hiattice, and to) is equal on all sites (see supporting in- 
formation for details). 

As probe light couples directly to the rotational motion 
of the dipoles, it is possible to measure the single spin-flip 
response of the system in order to detect and then charac- 
terize the FCI state. For example, the spectral function 
can be measured at finite energy and momentum using 
two-photon Bragg spectroscopy, providing direct infor- 
mation regarding fractionalization [30, 31]. On the edge, 
one should observe gapless chiral Luttinger liquid behav- 
ior, while in the bulk, the response should exhibit a gap 
to the multi-quasiparticle continuum. Such a gap mani- 
fests as an effective "magnetization" plateau as shown in 
Fig. 3D [11]. In this context, the fractionalization of the 
FCI state would lead to a nontrivial power-law onset of 
tunneling response. Alternatively, time- of- flight absorp- 
tion imaging can provide snapshots of the ground state 
liquid configurations. Such measurements allow direct 
exploration of the structure factor of the FCI wavefunc- 
tion and raise the tantalizing prospect of directly imag- 
ing fractionalization in real space. For instance, one can 
introduce a spin-flip into the bulk, allow the pair of frac- 
tionalized quasiparticles to propagate and finally, read- 
out using full density microscopy. 

In summary, we have introduced a scheme by which the 
collective physics of dipolar interacting spins can be har- 
nessed to realize a fractional Chern insulator. While we 
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have focused our discussion on polar molecules, our pro- 
posal can be realized in any system composed of electric 
or magnetic dipolar interacting generalized spins; such 
degrees of freedom are found in a diverse array of contexts 
ranging from magnetic atoms and Rydberg ensembles to 
solid-state spins [19, 32]. In particular, for exchange cou- 
pled electronic spin dimers or hyperfine coupled nuclear 
and electronic spins, one finds an effective level-structure 
nearly identical to that depicted in Fig. IB. The dipolar 
interaction between such coupled spins also yields topo- 
logically nontrivial, flat, spin-flip band-structures, en- 
abling the potential realization of a solid-state Chern in- 
sulator. While it is challenging to create electromagnetic 
gradients on the shortened length scales associated with 
such solids, it may nevertheless be possible to implement 
a flux lattice via homogeneous driving, staggered fields 
and spin-orbit interaction. With this in mind, several 
intriguing directions can be considered, including engi- 
neered solid-state defect arrays, magnetic atoms trapped 



in a quantum solid matrix, and dimer-based magnetic 
insulators [18, 19, 32]. Such approaches promise the ad- 
vantage of trapping dipoles at closer distances than in an 
optical lattice, and in the case of quantum solid matri- 
ces, such systems are already commonly used for preci- 
sion laser spectroscopy. With the enhanced interaction 
strength in a solid-state environment, strongly correlated 
many-body phases could be extremely robust. 
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I. DERIVING THE EFFECTIVE HAMILTONIAN 

Here, we derive the effective Hamiltonian, Hb = — Ylij ^ij a l a j + \ Vij n i n j (Eqn. 3 of the main text). The geometry 
is shown in the inset of Fig. 1 A. The molecules lie in the X-Y plane and the applied DC electric field has spherical coordinates 
(Go, $o) m this basis. To simplify the notation, we define |t) = s|l, — 1) + 1) + w|l, 0), where s = ^2^4/^, v — QiQs/Q,, 
w = Consider two dressed molecules at positions i and j separated by R = (R, 0, (j>) (spherical coordinates in the 

{x, y, z) basis). The dipolar interaction between these molecules can be written in spherical tensor form as: 

#dd = -4^^ ^(-l)^(^0)T g 2 (d«,d^), (SI) 

where (0, (j)) is the spherical harmonic of degree k and z angular momentum q in the normalization of [1]. Here, T 2 is 

the rank 2 spherical tensor generated from the dipole operators; in particular, 7± 2 (dW,dW) = d^±d±\ Tj^dW, dCO) = 

(d$d^ +dfd { £ ) } /v^,7?(dW,dW)) = (d^d^ +2d z i) d z j) H-d^d^) /V6,andd± = T(d x ±id y )/V2[l]. Expressing 

the dipolar interaction in this form allows us to isolate energy conserving terms. We assume that the energy difference between 
|1, 0) and |1, ±1) (A in Fig. IB) is larger than the scale of the dipole-dipole interactions. Under this assumption, T± x terms of 
the dipolar interaction are energy non-conserving and thus highly suppressed. 

We now consider the three resonant contributions to the hopping (fij) matrix elements [2-5], 

(Uj\n\UW = \j\[dl w* Wj -\dlMvi + s* Sj % (S2) 
{Uj\Tl 2 \UW = -d 2 01 (v* Sj ), (S3) 

<t4il^ 2 l^ti> = -<%Mvj), (S4) 

where doo = (1, 0|d^|0, 0) and doi = (l,±l|d±|0,0). Combined with the spherical harmonic coefficients of Eq. SI, these terms 
determine the directionally dependent hopping Uj of the spin flips. 

The interactions Vij between the spin flips arise as a consequence of the induced dipole moment which each molecule acquires 
in an applied electric field. Thus, Vij = (t4j| H dd |t4j) + (Uij\H dd \Uij) ~ (t4j| H dd |t4j) - {U\j\ H dd \U\j) can be 
calculated in a similar fashion [2]. First, let us define the induced dipolar moment of a molecule on site i, d^ i = d\(\si\ 2 + 
\vi\ 2 ) + no\wi\ 2 , where d\ = (1, ±l|d^|l, ±1) and /jLq = (1, 0|d^|l, 0) that of the |1, 0) state. The contributing terms to Vij are 
then (suppressing ij superscripts in d): 

(I4j|44 + ^(d+d_+d_d + )||4 j ) = d§, (S5) 

(tilj\d z d z + ^(d + d_ + d_d + ) |t4j) = d t .d , (S6) 
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(Utj\d z d z + -(d+d- + d-d+)\Utj) = d d tjJ (S7) 



(titj\d z d z + -(d + d_ + d_d + ) |t4j) = d ti d t . - -/^(s^i^s* + w^VjW* + c.c), (S8) 
(titj|d+d+ |titj> = -Moi(«t^>j^* +^<s^*), (S9) 
(titj|d_d_|titj> = -MoiK 5 *^^* +^>j5*), (S10) 

where do = (0,0|d 2 |0,0) and /iqi = (1, =bl|d±|l, 0) is the transition dipole moment between |1,0) and |1,±1). From H dd 
there also exists an on-site potential tu = Sj/;(44j| Il4j) — (t4j| Hdd |t4j)) which varies between sites; however, 
as we will see below, inhomogeneities in tu can be regulated using optical lattice tensor shifts. Finally, we note that we have 
dropped a uniform chemical potential term associated with the molecule's rotational constant 2B (see Fig. IB of the maintext). 
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To obtain topological flat bands, we adjust the optical radiation to generate four different types of sites {a, ft, A, B} (in relation 
to the M-scheme). By restricting the variation of |t) on (a vs. A) sites and on (b vs. B) sites, we ensure that both Uj and Vij 
are invariant under the direct lattice vectors #1, #2> enhancing the symmetry to that of a checkerboard lattice with a two-site 
translational unit cell. This small generalization from a two-site model [2] provides an important minus-sign freedom in the 
choice of w between lowercase and uppercase letter sites, which we exploit in tuning the Chern band structures. The freedom 
can be seen by examining the constraints imposed by requiring Uj and Vij to be invariant under a o A and b <-> B (i.e. 
translation by g2). In particular, the relevant constraints allow w a /b = w^/b or w a/b = ~ w a/b as solutions. 



II. EXPERIMENTAL IMPLEMENTATION IN 40 K 87 RB 



Molecular Hyperfine Structure: 



In this section, we consider the complications in our scheme due to the hyperfine structure of diatomic polar molecules such 
as 40 K 87 Rb [6-9]. The molecular rotational degree of freedom is naturally coupled to the nuclear spins, I\ = 4 and I2 = 3/2 
of potassium and rubidium. The hyperfine Hamiltonian is dominated by the nuclear quadrupole interaction, which has a typical 
strength Hq ~ 1MHz (for 40 K 87 Rb). This interaction splits the degeneracy between the |1, ±1) rotational states implying that 
our workhorse, the T± 2 terms of H^d, are off-resonant. To overcome this issue, one can simply ensure that the optical dressing 
ft (in the M-scheme) is much stronger than Hhf- This ensures that the hyperfine interaction is unable to couple the dark state to 
other dressed eigenstates. 

One final issue to consider is the particular choice of nuclear spin states. Since the composition of the dark state differs on the 
four types of lattice sites {a, ft, A, B}, molecules on these sites are subject to slightly different hyperfine potentials; in particular, 
(1a I H hf |t a) = (ta| H hf Ita) + (tB\H h f Its) = (t&l H h f |t&). Furthermore, the appropriate nuclear eigenstates will also 
depend on whether we are considering the ro vibrational ground state (|t)) or the dark state (|t)); this is because the decoupled 
nuclear spin basis (so-called Paschen-Bach regime [10]) is only valid in the first case. One can solve this issue by applying a 
static magnetic field ~ 10 3 G along the direction of the DC electric field, ensuring that the decoupled basis is appropriate for 
both \\) and |t) (note that such a field is already present in experiments such as [9]). Then, it only remains to choose a pair of 
nuclear eigenstates which have reasonable overlap and resonant energies. We have numerically verified that this is generically 
achievable. 



Optical Lattice Tensor Shifts: 



Similar to the hyperfine potential, A-type and 5-type sites feel different tensor shifts from the optical lattice. As alluded to in 
the main text, these tensor shifts can be exploited to compensate for dipolar induced tu terms. To start, let us consider a single 
optical field, E(R, t) = E(R)e~ luJt + h.c. which we use to create the lattice potential in the X direction. The optical potential 
is given by Hi at uce = —E(R)*a(uj)E(R), where E(R) = \E(R)\ J2 P Pp(R) e p, e p is the polarization basis, and a(uj) is the 
polarizabihty tensor of the molecule. Recasting the lattice potential in terms of spherical harmonics yields, 



Hj 



lattice 



-E 2 (R) 



2a± — ol\\ 



+ (ay -aa)^C^7 P 



(Sll) 



where a\\ is the polarizabihty along the internuclear axis, a± is the polarizabihty transverse to the internuclear axis, 70 

|A)| 2 - 1/3, 7±i = VV3 W± - /?;A)), and 7±2 = - VW^±- 
In our case, the optical lattice potential seen by ||) is, 



(i\ Hi att i ce |l) — —E 2 (R) 

while the potential seen by |t) is, 

(t| Hi a ttice |t) = 



2a± — a\\ 



+ (<*|| -cu)(0,0|C 2 |0,0) 7 o 



(S12) 



E 2 (R)[ 2a± + (ay - cu){7o(M 2 (1, -1| Cl |1, -1) 



+ M 2 (1, 1| Cl |1, 1) + M 2 (1,0| Cl |1,0» + 72^* (1, 1| C\ |1, -1) 
+ 7-2S*Wl,-l|C- 2 lM>}]- 



(S13) 



The energy difference SE = (f| Hi aUice It) — (II Hi attice \\) varies between A-type and B-type sites since the dressing pa- 
rameters {s, v, w} are site-dependent. The goal is to use this tunable tensor shift to compensate for dipolar induced tu terms. 
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FIG. SI: — Many-body Phases. (A) Depicts the spectrum associated with the knight's move solid phase in the reduced Brillouin zone (k x 
and k y are crystal momenta). The electric field tilt is Go = 0.05 and the field strength is \E\ ~ 32kV/cm. (B) depicts the spectrum associated 
with the striped phase (Go = 1.05 and \E\ ~ 28kV/cm). (C) depicts the spectrum associated with the striped supersolid phase (Go = 0.68 
and \E\ ~ 36kV/cm). (D) depicts the spectrum associated with the checkerboard phase (Go = 0.25 and \E\ ~ 40kV/cm). (E-G) depicts the 
structure factor of the KMS, CKB, SSS respectively for the same parameters as above (SS phase omitted). (H) depicts the spectral flow under 
magnetic flux insertion of the SSS phase. X-boundary condition twists produce little dispersion in the ground state energy while Y twists 
yield markedly superfluid behavior. Boundary conditions twists in both X and Y directions yield minimal superfluid response for the KMS, 
CKB and SS phases (not shown). 



Note that we can achieve propagation of the optical beams along any direction using only cr+ and tt light. Since our optical 
field never contains any a_ polarization, we find that j±2 terms are zero. Moreover, we have also dropped j±i terms, since 
A = E\$ — #1,1 ^> Hiattice- Combining the optical fields along the x, y and z direction, we have numerically verified that by 
simply adjusting the intensities of the lattice light, we can fully compensate for any inhomogeneous on-site potential. 

Finally, we demonstrate a simple configuration of Raman lasers (with wavelength Ao), which generates the M- scheme for 
the {a, 6, A, B} checkerboard lattice shown in Fig. 2A. We take the lattice constant to be A/, = Rq (Fig. 2A) and assume that 
Ao < Al; this can always be accomplished by increasing Xl (at the expense of weaker dipolar interactions). We can tilt the 
k- vectors propagating along X and along Y up or down out of the XY plane to give them a periodicity of Xl (in the XY plane). 
Similarly, we can tilt the fc-vectors propagating along (X ± Y) up or down out to give them a periodicity of V2Xl (in the XY 
plane). By using only four out of these eight beams and linearly polarizing them along k x z, we can obtain arbitrary Vt 2 and ^3 
on A and B sites with Qi = ^4 = 0. By utilizing the other four beams, we can obtain arbitrary Qi = SI4 on A and B sites (and 
their negatives on a and b sites) with Vt 2 — ^3 = 0. This immediately enables us to construct the four-site M-scheme. 



III. MANY-BODY PHASES 



Here, we provide a detailed description of the many-body phases which arise as one tunes the electric field parameters. First, 
we note that the field and dressing parameters for the phase diagram (Fig. 1C) are different than those for Figures 2B and 3. 
This is because the richest many-body phase diagram that we observe does not occur for the band structure with the largest 
flatness ratio. The band-structure depicted in Fig. 2B occurs at electric field tilt Oo = 0.68, $0 = 5.83, with optical dressing 
parameters: {0 aj #5, </> a , a a , a&, 7 a , 75} = {0.53, 0.97, 1.36, 3.49, 2.84, 2.03, 4.26, 3.84}, where we have parametrized: Si = 
sin(c^) sin(^), V{ = sin(a^) cos(^)e*^ and Wi = cos(c^)e* 7i . 

The phase diagram shown in Fig. 1C is computed by exact diagonalization for filling fraction v = 1/2 and a total of N s — 
24 sites (parameters: {9 , $0, a , 6 b , <j> a , </> h , a a , ot b , j a , lb} = {0.65,3.68,2.4,2.97,6.06,4.1,0.97,2.74,3.44,1.74}). The 
associated band- structure has an optimized flatness ratio, / « 7. The dressing parameters on A-type and B-type sites are 
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FIG. S2: — Variational Mean-field theory. Depicts the variational mean-field theory for K Rb as a function of applied electric field 
strength and tilt angle Go- We consider twelve variational ansatz's including all allowed solids up to a quadrupled unit cell. The superfluid 
ansatz is optimized with respect to its winding and a relative phase difference between A- and 5-type sites. 



identical to those on a and b sites with the exception that = tt + 7a and 7# = tt + 75 (ie. w a / b = —wa/b)- At weak 
electric fields, E < 8 kV/cm, diagonalization reveals the v = 1/2 Fractional Chern Insulator. By changing both the strength and 
tilt (0o ) of the DC field, one can map out a phase diagram containing both conventional and topological phases. To isolate the 
effects of long-range interactions, we ensure that at each tilt, the flatness ratio remains the same for all DC field strengths. This 
can be achieved by re-optimizing the dressing parameters for each field strength. Alternatively, this corresponds to ensuring that 
doo = doi as the field increases; experimentally, one can realize this by dressing the |1, ±1) states with a long-lived metastable 
excited state [6-9]. Numerically, we implement this constraint by taking S{ — > and V{ — » ^doo/^oi- 

There exist four crystalline phases at strong DC electric fields whose diagnostics we depict in Fig. SI. At low and intermediate 
DC field strengths, we observe a large superfluid region. This phase is characterized by a unique ground state (typically in either 
the (0, 0) or (— n, 0) momentum sectors) and a fluid-like real space structure factor. The homogeneity of the superfluid changes 
as one adjusts the DC field strength and tilt. At very weak fields \E\ < 4kV/cm, A-type and B-type sites are equally populated; 
however, as one increases the field strength, the anisotropy of the long-range dipolar interaction yields anisotropy in the structure 
factor. To verify the nature these non-topological phases in the thermodynamic limit, we perform a detailed variational mean- 
field study. As shown in Fig. S2, this study confirms not only the existence of these phases, but also the qualitative location 
of the phase boundaries. Moreover, in the KMS, CKB and SS solids, the mean-field energies match nearly identically with the 
exact diagonalization energies. This suggests that product state wavefunctions are valid approximations in these regimes and 
hence, that one can easily prepare finite density, low-temperature states. 
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